WIRELESS CAPACITY AND ADMISSION CONTROL IN COGNITIVE RADIO 



MAGNUS M. HALLDORSSON AND PRADIPTA MITRA 



Abstract. We give algorithms with constant-factor performance guarantees for several capacity 
and throughput problems in the SINR model. The algorithms are all based on a novel LP for- 
mulation for capacity problems. First, we give a new constant-factor approximation algorithm 
for selecting the maximum subset of links that can be scheduled simultaneously, under any non- 
decreasing and sublinear power assignment. For the case of uniform power, we extend this to 
the case of variable QoS requirements and link-dependent noise terms. Second, we approximate a 
' problem related to cognitive radio: find a maximum set of links that can be simultaneously sched- 

uled without affecting a given set of previously assigned links. Finally, we obtain constant-factor 
approximation of weighted capacity under linear power assignment. 
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1. Introduction 

How much communication can be active simultaneously in a given wireless network? This is a 
topic of major research effort. We address this question in a more generalized setting than previously 



q 

^ ■ considered, and give efficient algorithms that achieve good performance guarantees based on a novel 

mathematical programming formulation. 

In the capacity problem in wireless networks, we are given a set of communication links in 
^ . a metric space, each consisting of a sender-receiver pair, and we seek to find the largest subset 

I of links that can transmit simultaneously within the model of interference. We adopt the SINR 

■ model of interference where transmission over a link succeeds if the received signal at the receiver 

is sufficiently large, compared to ambient noise and interference from other transmissions. This 
model has emerged as a superior model for wireless interference patterns, as it is both analytically 
manageable, and reasonably realistic, especially in comparison to graph based models [Tl ] [26 t l29]. 
We assume that the powers have been pre-assigned to the links, based only on the length of the 
links. Having such simple assignments can be of great benefit in a distributed context. 

The basic capacity problem has been addressed in numerous recent works. Constant-factor 
approximation algorithms have been given for uniform power [S] and more generally for any non- 
^ ■ decreasing sub-linear power assignment (see Sect. 2 for definitions) |16j, and for arbitrary power 

5^ ! [H]- These results assume a uniformity of the links, both in their signal characteristics as well 

as their value. They also assume that no other wireless activity is affecting these transmissions. 
We aim to handle more general scenarios, allowing for heterogeneity in link characteristics and 
environment. In particular, we address three extensions: 

(1) (QoS) Each link has its own signal requirements and its own ambient noise term. 

(2) (Weights) Each link has an associated weight, and the objective is to maximize the total 
weight of the satisfied links. 

(3) (Admission control) Certain communication is already taking place, which cannot be inter- 
fered with (possibly for regulatory reasons). 

We discuss each of these extensions further. 

Quality-of-Service requirements. The SINR achieved at a particular link determines the data- 
rate achieved at this link, or the quality of service (QoS). Different links may have a minimum 
acceptable QoS requirements, for example if one link is used for video transmission and another for 



data transmission. In addition, the noise level at receivers may not be the same across the network. 
This practically motivated version of the capacity problem has not been handled in much of the 
previous research [51 [T5l [16j . We tackle this problem, both as an interesting problem in its own 
right, and as a stepping stone for the following problem. 

Cognitive radio and admission control. Given are two sets of links V and C. The goal is to 
find Q'^ C such that "P U Q can transmit simultaneously and the size of Q is maximized. We refer 
to this as the admission control problem. 

This problem naturally arises in at least two application areas. The first is the so-called "cognitive 
radio" , which has been the object of intense study of late (see [21 [331 121] and their many refer- 
ences). This area has gained great salience due to recent regulatory changes in wireless bandwidth 
management. Though the exact technological scenario for cognitive radios is still being figured out, 
the essential point is as follows: a wireless channel is allocated to a "primary user"; however, one 
would like to accommodate more users in the channel, as long as the primary user remains feasible. 
This clearly is an instance of the above mentioned problem, with V being typically small (in fact, 
perhaps, just 1). 

However, there is a more "classical" source of the same problem, referred to as admission control 
or access control [71 [36] , sometimes referred to as "active link protection" [5] . The capacity problem, 
is its basic form, captures a scenario where each slot is independent of previous slots. In practice 
however, links can require sustained communication (and different links for different periods of 
time). Thus, in certain applications, a more realistic model is to maximize capacity under the 
constraint that older links still communicating not be disturbed. This again, is exactly the problem 
defined above (but perhaps with a typically larger V). Though heuristic approaches to this problem 
abound, we are unaware of rigorous algorithmic results in the SINR model. 

Weighted capacity. In this problem each link v is associated with a non-negative weight and 
the goal is to find a feasible set OPT so as to maximize 'YIv&opt'^-v 

This weighted capacity problem is a natural extension of the capacity problem, and a case can 
be made for theoretical investigation for this reason alone. As it happens, though, the problem is 
further motivated by questions about stability in queuing theory. In this setting, packets arrive at 
network nodes according to some stochastic process, and the problem is to characterize the set of 
arrival rates under which the network can be stabilized, i.e., the network queues remain bounded. 
In the case of wireless networking stability, the seminal work of Tassiulas and Ephremides [M] 
established the existence policy that stabilizes the system under all arrival rates for which stability 
is potentially possible. This policy can be seen to be equivalent to solving the maximum weighted 
capacity problem in the SINR model. 

Solution method. It is easy to verify whether a given set of links is feasible. In fact, an appropriate 
power assignment that makes it feasible can be found efficiently. Namely, Eqn. [1] can be cast as 
a linear program with P^'s as variables, which can thus be solved optimally. Indeed, there is a 
large body of work where one starts with a feasible set and then tries to optimize over some other 
criteria, say to minimize the power consumed [5]. 

Naturally, one doesn't expect this approach to work directly for the capacity problem introduced 
before, which is "combinatorial" and in fact happens to be NP-hard [9]. What is perhaps more 
surprising is that the capacity problem does not appear to easily admit a linear programming 
relaxation either, even for simple cases. Most algorithms developed for the capacity problem have 
thus been very simple greedy algorithms [H [HI [21] , with some exceptions [171 S] ■ 

In this work, starting from a simple observation, we develop an integer program that approxi- 
mates the capacity problem for a large class of oblivious power assignments. We then show how 
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to round the corresponding linear programming relaxation to get a constant factor approxima- 
tion. Thus we recover the main result of [T6] but via linear programming as opposed to a greedy 
algorithm. We also show that the LP formulation can be easily modified to tackle a class of impor- 
tant problems where greedy algorithms do not appear to work very well, including the problems 
discussed above. 



2. Preliminaries and results 

The capacity problem in the SINR model is defined as follows. We are given a set L of n links, 
each consisting of a sender and receiver pair (^Sy^Vy^^ which are points in a metric space with a 
distance metric d. The asymmetric distance from link w to link v is the distance from w^s sender 
to v^s receiver, denoted d^v = d{sw,ry). Each link v has been assigned transmission power Py. A 
link V succeeds if 

(1) Py/dyy 

where is the ambient noise, /3 is the required SINR level, a > is the path loss constant, and 
S 3 V is the set of concurrent transmissions. A set S is feasible if the above constraint holds for all 
V £ S. Thus the capacity problem is equivalent to finding the feasible subset 5 C L of maximum 
size. 

Let £y = dyy denote the length of link v. Let A denote the ratio between the maximum and 
minimum length of a link. A power assignment P is non- decreasing if Py > Py, whenever > iy^ 
and sub-linear if ^ < whenever iy > iw We will restrict our attention to this class or particular 
assignments belonging to this class. Note that this class essentially contains all "natural" length 
based assignments, and specifically all well studied length based power assignments. These include 
the uniform power assignment, where all links use the same power, linear power assignment where 
Py = ly (which is thought to be energy efficient), and mean power assignment where Py = ly^"^ 
which is known to be essentially the "best" length-based assignment as far as capacity is concerned 

[laile]. 

Afjectance. We will use the notion of affectance, introduced in [HI [15] and refined in [22] to the 
thresholded form used here, which has a number of technical advantages. The affectance a^{v) on 
link V from another link w, with a given power assignment P, is the interference oi w on v relative 
to the power received, or 

a^{v) = mm< l,Cy— ■ — 

where Cy = /3/(l — f3Nyi" / Py) is a constant depending only on the parameters of the link v. 

We will drop P when clear from context. Let ay{v) = 0. For a set S of links and a link v, let 

o-viS) = E«;e5"f(^) (^s{v) = Et«65"«'(^)- For sets S and R, aR{S) = YlyenT^ueS ^^^i^)- 
Using such notation, Eqn. [T]can be rewritten as follows, which we will adopt: 

(2) asiv) < 1 

In the variable QoS version of the capacity problem, /3 and A^ are no longer constants, but can 
be different for different links. Note that the definition of affectance stays the same apart from 
a changed definition of Cy = /3y/{l — l3yNy(.'^ / Py) where fiy and A'„ are respectively the signal 
requirement and noise level for v. 

For all problems that we consider, we will OPT to mean the optimal solution, which will apply 
to the problem being discussed at that point. 

3 



Our results. We prove the following results. 

Theorem 1. For length monotone, sub-linear power assignments, there is constant- approximation 
algorithm for the wireless capacity problem. For uniform power, there is a constant- approximation 
algorithm in the QoS generalization. 

The first part (not involving QoS) is the same as the main result proven in [16J, but via a linear 
programming relaxation. 

Theorem 2. For the admission control problem with uniform power, 

a) There is a 0{\V\) approximation algorithm. 

b) If the optimum solution \OPT\ > 7i|'P|\/log {Vl (for some constant there is a constant- 
approximation algorithm. 

Specifically, for the "cognitive radio" case of the problem where V = 1 (or small, at any rate) 
we get a constant factor approximation for uniform power. There is no straight-forward greedy 
algorithm to tackle this problem. We believe that a greedy algorithm for the variable QoS problem 
is possible, but even if this is true, the resultant version for admission control would result in 
approximation factor worse than our results by a O(logn) factor. Additionally, we see no way of 
utilizing the OPT > |7^| y^log \V\ condition in the greedy algorithm. 

Finally, 

Theorem 3. For linear power, there is constant- approximation algorithm for weighted capacity 
problem. 

For this problem, greedy algorithms combined with some basic observations can yield a 0(min{log A, log n})- 
approximation algorithm (we describe this algorithm in detail in Section [6] when we experimentally 
compare it with our LP based algorithm). 

We remark that our results holds in arbitrary metric space, independent of the path loss constant 
a, and faithfully treat the ambient noise term. 

Related Work. The first work to study capacity of randomly deployed networks was the work 
of Gupta and Kumar |12] . Rigorous worst case algorithmic analysis started with the work of 
Moscibroda and Wattenhofer [27] . who studied of the scheduling complexity of arbitrary set of 
wireless links. Early work on approximation algorithms produced approximation factors that grew 
with structural properties of the network |30l [281 E] • 

The first constant factor approximation algorithm was obtained for capacity problem for uniform 
power in [8] (see also [15]) in with a > 2. Fanghanel, Kesselheim and Vocking [6] gave an 
algorithm that uses at most 0{OPT + log^n) slots for the scheduling problem with linear power 
assignment P^ = that holds in general distance metrics. 

Recently, Kesselheim obtained a constant-approximation algorithm for the capacity problem 
with power control for doubling metrics and O(logn) for general metrics [21]. In another work |16j . 
constant factor approximation was achieved for all non-decreasing, sub-linear power assignments. 
The greedy algorithms of [HJ [T5l [T6] can be modified to handle the problems we address here, and 
these algorithms essentially constitute the previous best results on these problems. 

As far as we can ascertain, the algorithmic situation for the admission control and weighted 
capacity is somewhat similar to the situation the "basic" capacity problem was in before the array 
of results mentioned above. Thus, we have a large body of works and results in different settings 
motivating the study of these questions, but no worst case algorithmic results. 

The works on the emergent field of cognitive radio are too numerous to adequately cover. We 
refer the reader to [2] for a thorough discussion. For results on capacity of networks in a cognitive 
radio context see [191132] etc. ("capacity" not necessarily meaning the exact same thing we do). For 
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the stability problem in a queuing theory setting that gives rise to the weighted capacity problem 
there are many works in graph based models |31l [20l [25] as well as recent ones on the SINR model 
[23] • The weighted capacity problem was recently studied in ^35j, where the authors propose a 
version of the greedy-based O(logn) approximation. 

In terms of using a LP approach in the SINR setting, there is recent work of Hoefer et. al. 
[17] . using related insights in their formulation. In the context of throughput maximization, [3] 
employed a linear programming solution. Being based on unit disc graphs, that approach does not 
lead to performance bounds we seek here. 

3. The basic capacity problem 

Let us first consider how one would attempt to write an Integer program (and a subsequent 
Linear programming relaxation) for the capacity problem. If the variable 5^ G {0, 1} denotes that 
link V was selected in the solution, we see that for selected links, the condition YlueL '^uiv)5y < 1 
would have to hold. This is quite nice and linear, except for the fact that we would have to somehow 
indicate that this condition need only hold for {v : 5v = 1}, and that no condition need hold for 
links in {v : 6^ = 0} . There appears no way to do this in a linear program. 

For clarity, we will first present our linear program (and the whole algorithm) below, and then 
in proving its correctness, we will describe how out algorithm evades the problem elucidated above. 

3.1. Algorithm. Our algorithm has three main steps: 

• Linear Program The first step of our algorithm is to solve the following linear program, 
with variables 5u, one corresponding to every link i^- Let C be a large enough constant. 

(LP) maximize (5m subject to 

(3) av{u)5y < C; Vn 

v/v>eu 

(4) Yl (^u{v)5,<C- Vn 

v,ev>iu 

(5) < 5„ < 1; Vn 

• Rounding We then "round" the fractional solution to this linear program in two steps. 

Let LP* be the value of the (fractional) solution to LP. 

First, we select a set i? = {n : = 1}, defined by binary variables s^, which are generated 
independently at random such that P(su = 1) = 5u (and thus P{su = 0) = 1 — (5„). 

Next, we choose a subset of R named S defined by = {n : = 1}, where s'„ is a binary 
random variable corresponding to this second round of selection. The variable s!^ is defined 
as follows: = 1 iff = 1 and the following two conditions hold: 

(6) ay{u)sv < 3C 

v,£v>eu 

(7) Yl ('uiv)s,<3C 

• Final Selection Finally, a feasible set is extracted from S using a simple signal- strengthening 
technique which we will detail later. 

3.2. Analysis. We need the following definitions. 

Definition 4. A link set L is 7-feasible (resp., 7-anti-feasiblej, if aiiu) < 7 for all u £ L (resp. 
if au{L) < 7 for all u £ L). A link set is 7-bi-feasible if it is both ^-feasible and ^ -anti- feasible. 
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We will simply write "feasible", "anti- feasible" and "bi- feasible" when 7 = 1. 
Our first step is to show that the solution to the linear program is an approximation to the 
capacity problem, or more formally: 

Lemma 5. Let LP* be the value of the optimal solution of LP. Then, LP* = ^}{\OPT\). 

Proof. To prove this, it suffices to construct a solution 6u (for all u) to the LP such that ^„ 6u = 
0,{\OPT\), and that satisfies all the constraints in the linear program. 

Since OPT is feasible, there is a 2-bi-feasible subset W such that \W\ > \OPT\/2 (See [13] for 
a simple proof of this fact). 

Now construct the solution by setting 6u = ^ u £ W and 6u = otherwise. Thus J2u ~ 
\W\ > \0PT\/2. Lemma[5]can be proven then if we can show that Conditions [3] and H] hold for this 
solution, and thus form a valid solution to LP. These follow directly from two Lemmas noted below 
(Lemmas [6] and [71) , by setting C to be larger than the implicit constants in those two Lemmas. □ 

The following Lemma was proven in [22]. For completeness, we give a proof in the appendix that 
holds for arbitrary ambient noise. 

Lemma 6. Assume L' is ^-feasible using a non- decreasing, sub-linear power assignment. Let u be 
any link such that 1^ < H-v for all v G L' . Then 0^1 (u) = 0(7). 

The next Lemma, something of a dual of the previous one, was proven recently in [14j: 

Lemma 7. Assume L' is ^- anti- feasible using a non- decreasing, sub-linear power assignment. Let 
u be any link such that < iv for all v S L' . Then au{L') = 0(7). 

Remarks: Lemmas [6] and [7] hold the crucial insight that allow us to circumvent the problem 
mentioned at the beginning of this section. Note how these lemmas bound the affectance to and 
from a link u without the condition that be a part of the feasible (or anti-feasible) set L. This 
allows us to evade the issue of having to express conditions that only apply for links in the solution 
set. Instead we can write constraints (Equations [3] and [3|) which apply to all links. 

The next step is to analyze the Rounding phase. In particular, we claim that 

Lemma 8. ^{\S\) = n{LP*) = n{\OPT\). 

Proof. Recall that 5 = {u : = 1}. Then by linearity of expectation. 



E(|5|) = E{J2s'J = J2n<) 

u u 

= J]P« = 1) = ^P« = 1|S„ = 1)P(S„ = 1) 

u u 

(8) = j;P« = l|5„ = l)5„ . 

u 

where we use P(s„ = 1) = 

Let p{u) denote the indicator random variable of the event that both Cond. [6]and[7]are fulfilled 
for link u. Then P(s^ = = 1) = P(yo(u)|s„ = 1). The point to note here is that the events p{u) 
and = 1 are independent since the random variable Su is not involved in the former (because 
au{u) = 0). 

We win prove below that P(su = 1) >\ (Lemma [9]). 

Thus, P(s; = = 1) = P« = l)>\. Now continuing with Eqn. [HI E(|5|) > \ Y.u^u > 
\LP*. 

□ 
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As promised, we lower bound P(s^ = 1): 
Lemma 9. P« = 1) > | 

Proof. By Eqn. d it holds that ^{Yle,>e^ ay{u)sv) = T,e,>e^ay{u)E{sy) = Y^i^>i^ay{u)6y < C. 
Thus by Markov inequality, the probability that Cond. [6]fails is at most 1/3. The same applies to 
Cond. [3 using Eqn. HI The Lemma then follows by the union bound. □ 

Finally, we need to show that we can extract a large feasible subset from S in the Final Selection 
phase. The following signal- strengthening lemma from [15j will be frequently useful. 

Lemma 10. [Thm. 1 of |15j . slightly restated] If is an ^-feasible set, then S-y can partitioned 
in to 6-feasible sets, for any 6 < j. 

Lemma 11. There is an efficient algorithm to find a feasible set S" C such that \S'\ = il.{\S\). 
Proof. By conditions [6] and [71 the average affectance in the selected set S is 

^ a^,(n) = ^ ^ {a^{u) + au{v)) 

Define Si = {u £ S : as{u) < 12}. From the above bound on average affectance, it is easy to see 
that l^il > |5'|/2. Finally, by the signal strengthening lemma, we can find a feasible set 5*2 such 
that l^sl = = □ 

The first part of Thm. [1] now clearly follows. The part of Thm. [1] about uniform power in the 
variable QoS case will be handled in the next section. 

The algorithms in the following two sections will follow the same tri-partite design of LP, Round- 
ing and Final Selection. Due to space constraints, we will mostly focus on the changes in the LP 
formulation, and when appropriate, the changes in the Rounding phase, without proving everything 
from scratch. 

4. Cognitive radio/ Admission control 

Variable noise and signal requirements (QoS). Recall that in this variation of the problem 
each link v has a separate QoS 13^ and noise level and definition of affectance changes accordingly. 
If a link set is such that Py < ciPu for all u ^ v for some unspecified constant ci, we call the power 
assignment nearly uniform. The following holds. 

Lemma 12. Assume L' is anti-feasible and u is some link. Assume that all links use a nearly 
uniform power assignment. Then au{L') = 0(1). 

The proof is a standard modification of the same result for uniform power with constant f3 and 
(see, for example. Lemma 11 of ^j). Our proof of Lemma [6] provided in the appendix gives a 

general idea of this type of proof, and we mention after that proof the main changes needed to 

achieve Lemma [T2l 

The following modified LP can be used for uniform power capacity in this setting: 

(LP2) maximize 5u subject to 

(9) ^ 

< 5n < 1; V-u 
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The additional steps after solving the LP and the analysis follows the same lines as Thm. [T] 
(which we omit due to space constraints). 

Admission control. Now we can focus on the admission control problem for which we will use 
some ideas from the variable QoS case. 

We will prove the following more general result first. 

Theorem 13. Assume links in C use a nearly uniform power assignment. Assume that links in V 
use some arbitrary power assignment. Then we can approximate the admission control problem up 
to a factor of 0{\V\). 

Proof. Recall that the goal is to find OPT C £ of maximum size such that OPT U 7^ is feasible. 
Thus in choosing OPT, we have to be careful about the affectance of V on OPT and vice-versa. 
Our approach is to handle the affectance from V as noise. In this regime, the new "noise" present 
at each link is the original noise N , plus the interference received from all links in V. Specifically, 
for M G U £, we define a variable noise level Nu = N + YliveV u^v Define a„(w) to be the 
affectance taking this variable noise into account. 
Now consider the following LP relaxation: 

(LP2) maximize 5u subject to 

(10) J] J]a,(w;)5, < |P|; 

(11) ^au{v)5^<C] \/ueC 

< 5n < 1; Vn G £ 
We show that the solution of LP2 is close to OPT. 

Lemma 14. Let LP2* be the value of the solution to LP2. Then LP2* = n{\OPT\). 

Proof. Consider a 2-anti-feasible subset of OPT, call this O. Consider the following solution to 
the LP: set (5^ = 1 if « G O and otherwise. Cond. [10] is satisfied since the incoming affectance 
on each link in V from OPT (and thus from O) is at most 1. The case of Cond. [TT] follows from 
Lemmaini Thus LP2* > \0\ = n{\OPT\). □ 

The next step is to round the fractional solution achieved from solving the LP. As before, we 
first set s^j = 1 with independent probability Let us define the event A as the condition 
SweP Siie£ ^''(^)^f — holding. Let us define the event Bu, for each link n G L, as the 
condition J2vec ^u{v)sv < 4C holding. 

We derive another round of selections by setting = 1 iff = 1 and both Bu and A occur. 
Thus, 

E(4) = P« = 1) = P(s„ = l^Bu^A)= FiBu AA\su = l)Su 
Now F{Bu A A\su = 1) _> 1 - F{Bu\su = 1) - ^A\su = 1). As we have seen before, Bu is 

independent of Su, thus ¥{Bu\su = 1) < 5 (via Cond. [TT] and Markov's inequality). 

On the other hand, A is not independent of Su- However, A occurring given Su — 1 is the 

same as Y.v(^c,vy^uYlw&r o.y{w)sy < 5|7^| - Y.w&s^y-i'^) being true. But Et«eT'0«M < 1^1 > by 

the definition of affectance. Thus P(^|s„ = 1) < 'F{J2v(^c,v^uJ2weV^^(''^)^v > 91"^!) ^ j- Thus 

finally, F{Bu A A\su = 1)>1-|-|> 0.5. Therefore, E{s'J > 0.5(5„. 
After the last round of selection, we thus get a set C such that 

• R = 0,{LP2*), in expectation 
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• J^wf^R^viw) < 4C for all V e R 

Using averaging arguments and signal strengthening as before, we can extract R' = Q{R) which is 
feasible. To complete the solution, we need to extract a subset of R' such that 

(12) ay{w) < 1 for all w £ V 

From the condition X^^gp ai;(w)s't, < 10|7^|, it is not hard to see that the set of selected 

links from L can be partitioned into 0(|P|) sets such that Eqn. [12] holds. This gives us the sought 
after 0(17^1 )-approximation. □ 

Thm. [13] implies part a) of Thm.[2]directly. We note that this implies a 0(1^1 log A)-approximation 
algorithm that holds under any other non-decreasing sublinear power assignment, by partitioning 
the linkset into at most log A sets of nearly-uniform power. 

We prove the last part of Thm. [2] below. 



Theorem 15. Let k = \V\. If \OPT\ > jik^/logk, for a large enough constant 71, then there is a 
constant-factor approximation algorithm for the admission control problem for uniform power. 

First we show that if OPT is large, we can assume that the affectances from OPT to V are 
small. 

Lemma 16. Assume OPT > 7iA;\/log k, for a large enough constant 71. Define L' = {u G C : 
a^{v) < j^^^^ for all v £ V} and OPT' = L' n OPT. Then OPT' = n{OPT). 

Proof. To see this, note that aopT(P) < k, since V must be feasible in presence of OPT. Now, 
defining OL = OPT\OPT', aoUV) > \OL\j^^^. Thus \OL\j^^ < k, or \OL\ < lO^hik-k 

and finally \OPT'\ > \OPT\ - \OL\ > (71 - 10) ^/logl • k = n{\OPT\) if 71 is large enough. □ 

We can also claim a strengthening property. 

Lemma 17. Assume R is a set such that for all u £ R, Ouiv) < I , for all v G V, and 

o-r{v) < 1 for all v £ V. Then there is a subset R with \R\ = Q{\R\), such that a^(£) < 1/3 for all 
i £ V and such a subset can be found in polynomial time with high probability. 

Proof. Simply select each link in R with probability ^. Let the set of selected links be Q. Then 
E(aQ(u)) < I for all v £V. Consider a fixed v £ V. Now aQ{v) = ^ugr'^u{v)su where Su is the 
iid random variable indicating selection in to Q- We can use Hoeff ding's inequality to get a large 
deviation bound. 

Theorem 18. [Hoeffding, [IS]/ Let the independent random variables Xi, . . .Xn be bounded, i.e., 
Xu £ [bu,du], and let X = Yj^Xu. Then, 

2t^n'^ 

Set Xu to be au{v)su- We can verify that given our assumptions, setting 6„ = and du — ^ 
suffices. Setting t = h, 



'{X - E{X) >t)< exp 



lOVlosfc 



aqiv) > ^) < F [agiv) - E{aQ{v)) > ^ 



This imphes, by the union bound, that with probabihty at least aQ{v) < | for all v V 
simultaneously. We now have proof of not only the existential statement, but the algorithmic one, 
since we can repeat the random experiments multiple times to get the high-probability result. □ 

Note that the above holds equally for any affectance function (specifically, the case a). 

Now we can describe the linear programming relaxation. First note that by virtue of Lemma 
[TBI it suffices to assume the input instance is L' and the optimum is OPT' (links not in L' can 
be thrown out by simple pre-processing). Let us reuse notation C and OPT to refer to this new 
instance after pre-processing. 

Consider the following linear program {LP4) 

(LP4) maximize subject to 

(13) ^viw)5^ < 1/3; Vu; G P 

(14) Yau{v)6.a<C; £ C 

vec 

0<Su<l; yu£ c 
We claim that this is a relaxation up to constant factors. 
Lemma 19. Let LP* be the optimal value of LP4. The LP* = n{\OPT\). 

Proof. (Sketch) By Lemma [T71 there exists O C OPT such that aoiw) < ^ for all w ^ V. Now 
consider a 2-anti- feasible subset O' of O. Consider as the solution of the LP, (5^ = 1 if n G O' and 
6u = otherwise. That Cond. [13] is satisfied follows from the claim that ao{w) < \ for all w ^ V. 
Cond. [T3] follows from anti- feasibility arguments along the lines made before. □ 

We can round this solution in the same way as before, with a two stage selection process. The 
proof varies only in that we need to claim that after the second selection (characterized by Bernoulli 
variable s^), Ylv&c'^vi^)^v — ^ with high probability, simultaneously for all w gV. This follows 
from an argument similar to Lemma [T71 using the fact that affectances are bounded by iQ^^ogk ^^'^ 
using the Hoeffding's inequality. 

5. Weighted capacity 

Recall that the result for weighted capacity applies only to linear power. For linear power the 
following stronger version of Lemma [6] holds. 

Lemma 20. Assume L' is feasible using linear power and £ is any link (also using linear power). 
Then, aL'{i) = 0{l). 

The proof is nearly identical to that of Lemma [6l as elaborated in the appendix. 
We can now write the following LP relaxation for the weighted capacity problem. 

(LP5) maximize WuSu subject to 

(15) ^^at,(n)(5^ < C; Wu 

< 5„ < 1; Vn G L 

Proof, of Thm. [3] (sketch) The proof is rather like that of Thm. [TJ Once again, we select each link 
into a set R with probability 6u (characterized by Bernoulli variable Su for each link u) and then do 
a further selection by setting S = {u : s'^ = 1} where s'„ = 1 iff = 1 and J2veR'^'"^'^) — 
in the proof of Thm.[Tl one can show that E(s'„) = 0,{6u) and thus the expected weighted output is 
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^(X^u '^u^u) which is within a constant factor of the optimum. Finahy, the set 5* can be partitioned 
into a constant number of feasible subsets using signal strengthening (Lemma llOp . completing the 
proof. □ 

For other power assignments, such as uniform power, we seem to be within striking distance of 
a 0(l)-approximation. This is unfortunately not the case, we can only claim a poly-logarithmic 
approximation, worse than the greedy case. However, as we show in Section [6l in practice the LP 
approach might be applicable to these other power assignments as well. 

6. Simulations 

In this section, we present results from simulation experiments. We focus on the weighted 
capacity problem for our experiments. It is difficult to conduct a comparative experiment for the 
admission control problem, there being no obvious previous algorithm to compare it with. 

In contrast, the weighted capacity problem admits straightforward modifications of to the greedy 
algorithm, and thus a better comparative benchmark for our algorithm. Two natural greedy algo- 
rithms can be proposed: 

• Using weight classes: Let max^ u;^ = n (by scaling). Now we can assume that all 
Wu £ This is because links with smaller weights can be discarded without losing 
more than a factor of 2 in the approximation quality. Now divide the links into log n weight 
classes, the weig ht class Wt is defined by Wt = {u : Wu £ [2*, 2*+^]} for t = to log n - 1. 
Now if we consider links belonging to a single Wt, the weights do not matter (up to a factor 
of 2). We simply run the greedy algorithm of |16j for each Wt, and output the solution for 
the best weight class. This gives a straightforward O(logn) approximation factor. 

• Using length classes: Let min^^^ = 1 by scaling and let max-^ — Divide the 
links into length classes Lt = {u : £u ^ [2*, 2*"^^]} for t = to log A. Within Lj we can 
choose to run the greedy algorithm on the links in any order since the lengths are essentially 
the same, thus we go through links according to descending order of weights, achieving a 
constant factor approximation on Lt- We choose the solution for the best Lt, thus getting 
a O(logA) approximation. 

Thus, comparing the two greedy algorithms, we achieve a 0(min(log n, log A)) approximation. In 
what follows, we shall refer to this joint algorithm as "greedy algorithm". 

Experimental setup. We randomly generated the instances. Some important parameters of the 
experiments are as follows: 

(1) A: The maximum length of a link (the implicit minimum being 1) 

(2) R: A number indicating that the sender of a link is chosen from a R x R square 

(3) n: number of links 

We also use N = 0, /3 = 1, and a = 2.5. The instances were generated as follows. For each link, the 
sender was chosen randomly from a Rx R square. The length of the link was chosen randomly from 
[1,A]. The receiver was thus placed at this distance from the sender and at a random direction. 
The weight was chosen independently from [l,n]. We will mention different weight distributions 
later, and mention this initial choice of weight distribution as the ordinary distribution. 

One crucial aspect of both greedy algorithms as well as the LP algorithm is the constants used. 
For the LP algorithm, this is the constant C in Eqn. [T5j The greedy algorithm of [16] also depends 
on a constant. Though theoretical bounds for these constants are available, it has been observed 
before that these theoretical bounds do not perform the best in practice [1]. We run all algorithms 
with different values of the constant in question, running over reasonable values in small increments, 
and choosing the best solution for each algorithm separately. We ran our experiments in MATLAB, 
and used the convex optimization package CVX |10j to solve the LP. 
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500 1000 2000 5000 

R 

Figure 1. Simulation with linear power and n = 400 (ordinary weight assignment). 
Individual lines refer to different values of A. The ratio of the solution from the LP 
algorithm to the greedy algorithm is plotted on the Y-axis against the density. 




Figure 2. Simulation with linear power and n = 600 (ordinary weight assignment). 

The overall message from the experiments is that using the linear programming formulation gives 
a substantial improvement in the solution quality in many cases. On the other hand, the greedy 
algorithm is also not without merit, and can outperform the LP in certain other situations. It 
appears that the smaller the maximum feasible set is, the better greedy does, while as the solution 
size/quality improves, LP outperforms greedy. This is not surprising. When the set is really dense 
and the link lengths are large, the quality of the solution is bad and the cost incurred by greedy 
due to length-class or weight-class partitions is minimal. 

In Fig. [H we see the results for linear power with 400 links. On Y-axis is plotted ^j-, where 
SLP and SG are, respectively, the quality of the solution found from the linear programming 
algorithm and the greedy algorithm. As alluded before, the greedy algorithm does better when A 
and density are both large (these are the points for which the Y-axis value is lesser than 1), with 
the trend reversing when these change. Running the same experiment run for an increased number 
of links n = 600 confirms these trends (Fig. [2]) . 

We experimented with different distributions on the weights, to see if changes here change the 
solution trend significantly. We tried the following weight distributions. 

• Reversed: Set the weight of link f to be 1 /w^ where is chosen according to the ordinary 
distribution. 

• Length determined: Set weight of the link to be equal to its length. 
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Figure 3. Simulation with linear power and n = 600 and Reversed weight assignment. 




Figure 4. Simulation with uniform power and n = 400 and the ordinary weight assignment. 



• Weight class: Choose a parameter t randomly from [l,logn] and set weight to 2*. 

The overall trend is similar. For reversed and length determined distributions, LP did extremely 
well, whereas for the case of weight class distribution, greedy did much better, with LP only barely 
outperforming it in a few cases. This further points to the benefit of combining these algorithms 
in practice. The results for the reversed case are shown in Fig. [3j 

Next we experimented with uniform power. As we discussed in Section [5l for uniform power 
we can only claim a poly-logarithmic approximation factor. However, the bounds are only so bad 
on rather pathological instances and one needs to do some work to come up with them. Thus in 
practice, it is reasonable to assume that an LP approach will be not without benefit. This is indeed 
borne out by our experiments, as seen in Fig. SI 
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We give a proof of Lemma [6l originally due [22] , that holds also in the presence of arbitrary 
noise. 

Lemma El' If L is 7- feasible using a non-decreasing sublinear power assignment and u is a link 
such that iu < for all v £ L, then ai^u) = 0(7). 

Proof. Assume that L is a l/3"-feasible set. By the signal strengthening (Lemma llOp . this affects 
only the constant factor. 

Consider the link v G L such that d{ru,ry) is minimum. Also consider the link w £ L with 
d{sw,ru) minimum. Let D = d{ru,ry). We claim that for all links x m L, x ^ w, 

(16) d{s^,ru) >]^D . 

To prove this, assume, for contradiction, that d{sx,ru) < \D. Then, d{sw,ru) < ^D, by definition 
of w. Now, again by the definition of v, d{rx,ru) > D and d{rw,ru) > D. Thus 1^ > d{ru,ryj) — 
d{sw,ru) > y and similarly > y. On the other hand d(s^, s^,.) < d{su,,ru) + d{sx,ru) < "j + y < 
D. Now, ■ dxw < (4 + d{sw,s^)){iuj + d{syj,Sx)) < (4 + -D)(4 + ^) < 9^«)4, contradicting 
the following: 

Lemma 21 ([13]). Let u,v be links in a 1 / q"' -feasible set. Then, duv • d^u ^ ■ 44- 

Consider now any link x m L, x ^ w. By the triangle inquality and Eqn. 1161 dxv = d{sx,ry) < 
d{rv,ru) + d{sx,ru) = D + d{sx,ru) < 3d{sx,ru) = M^u- Now a^iu) < Cu-^^- Since 4 < 4, it 
holds that Cu < and by sub-linearity it holds that Pu/^u — Pv/^v- Thus, 

p pa "V^P 

(17) ax{u) < Cv-^-^ < = 3 a^ix), 

where the final equality follows from the feasibility of L. Finally, summing over all links in L 

aUu) = ^ ax{u) = a^{u) + ^ ax{u) 
xeL x<^L\{w} 

<l + 3° a«W < 1 + 3" -7 = 0(1) , 

x^L\{w} 

since Yl,xeL\{w} ^x{u) < aL{u) < 7 by assumption. □ 

Proofs of other lemmas: When using linear power, it holds for all links u and v that Cu = and 
the signal received Pu/^u — Pvl^'i-, satisfying Eqn. [T7] holds without the condition 4 < 4- This 
yields Lemma [20l 

By inverting the role of senders and receivers, we can obtain similar bounds on out-affectance 
(a„(L)) as above on in-affectance (aL(n)). IVIodulo this change, the proof of Lemma [7] is nearly 
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identical to the above proof of Lemma [6l and the argument for Lemma [T2] mirrors that of Lemma 
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